Melting behavior of ultrathin titanium nanowires 
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The thermal stabihty and melting behavior of ultrathin ti- 
tanium nanowires with multi-shell cylindrical structures are 
studied using molecular dynamic simulation. The melting 
temperatures of titanium nanowires show remarkable depen- 
dence on wire sizes and structures. For the nanowire thin- 
ner than 1.2 nm, there is no clear characteristic of first-order 
phase transition during the melting, implying a coexistence of 
solid and liquid phases due to finite size effect. An interesting 
structural transformation from helical multi-shell cylindrical 
to bulk-like rectangular is observed in the melting process of 
a thicker hexagonal nanowire with 1.7 nm diameter. 

61.46.+W, 68.65.+g, 73.61.-r 

Metal nanowires are current focus of intensive research 
due to both the fundamental interests of low-dimensional 
physics and the potential applications in nanoscale ma- 
terials and devices. Experimentally, stable ultrathin 
metal nanowires with diameter down to several nanome- 
ters and suffLeient length have been fabricated by var- 
ious mcthodscTEl. For example, Ti nanowires of a few 
nm in width were produced by Ar+ ions irradiation ori 
Ti thin layer by using carbon nanotube as a masko. 
Takayanagi's group, has successfully produced suspended 
stable gold wiresBB. Novel helical multi-shell structures 
were observed in those ultrathin gold nanowiresB. As 
compared with bulk materials and atomic clusters, one- 
dimensional (1-D) metal nanowire have both the 1-D pe- 
riodicity along wire axis similar to bulk solids and the 
larger surface to volume ratio like clusters. Therefore, 
the metal nanowires are expect to exhibit some unique 
features different from either bulk solids or nanoclusters. 

On theoretical side, the structures and properties of 
free-standing ultrathin petal nanowires have been stud- 
ied by several groupstjilZl. Based on genetic algorithm 
simulations with empirical potentials, we have systematp 
icalLy studied the helical multi-shell structures in AuEJ, 
TiEj, Zrll3 nanowires. The corresponding vibrational, 
electronic, and transport properties were also discussed. 
However, little is known about the thermal stability and 
melting behavior of the 1-D metal nanowires, especially 
for those transition metal nanowires with novel cylindri- 
cal multi-shell structures. Thus, it is important to clarify 
the relation between the melting behavior and the size 
or geometry ofrthese ultrathin metal nanowires. The fi- 
nite size effectst3~E3, e.g., solid- liquid coexistence, on the 
phase transition of such quasi-lD systems are also in- 
teresting. For thicker crystalline lead nanowires, surface 
melting phenomena associated with formation of a thin 



skin of highly mobile surface atoms were observed^. More 
recently, the multi-shell gold nanowires and their melt- 
ing behaviors have been studied by MD simulationsEj. 
In our previous work, the atomic structures of titanium 
nanowires with diameter from 0.7 to 1.7 nm have been 
optimized and helical multi- walled cylindrical structures 
have been obtainedEj. In this paper, we report molecular 
dynamic simulations on the thermodynamic properties of 
the titanium nanowires. The effect of nanowire size and 
geometry on the melting behavior will be discussed. 

In our simulations, Ti nanowires of sufficient length 
are modeled by supercell with 1-D periodical boundary 
condition along the wire axis direction. As a reason- 
able compromise between keeping the helicity in z axis 
direction and avoiding the nanowire breaking into clus- 
ters upon relaxation, the length of the supercell is cho- 
sen to be 1.256 nm. The interaction between titanium 
atoms is descrLbed by a well-fitted tight-binding many 
body potentials. To simulate the melting behavior of 
nanowires, we employ the constant temperature molec- 
ular dynamics (MD) method by HoovercJ, pticii have 
been extensively used in our previous worksE3 c3. The 
MD time step is chosen as 2.15 fs. At each temperature, 
the initial 10^ MD steps are used to bring the system 
into equilibrium. Then we monitor the internal energy 
E, root-mean-square (rms) fiuctuation of the interatomic 
distances (5, and heat capacity from the thermal statis- 
tical averages in the equilibrium canonical ensemble. At 
each temperature, 10^ MD steps are performed to record 
the thermodynamic average of these physical properties. 
The constant-temperature MD simulations start from a 
low temperature (400 K). The temperature gradually in- 
creases towards high temperature (1400 K) by 50 K per 
step. 

We first pjpr form comprehensive genetic algorithm 
simulationsliJIlj on the equilibrium structures of tita- 
nium nanowires. Several representative structures are 
presented in Fig.l. These stable multi-shell structures 
are composed by coaxial cylindrical shells, which have 
been theO|i;etiGall*-|predicted for Al, Pb, Zr, Cu, and Au 
nanowires ^ I13l3'E3 and experimentalhj— observed in Au 
nanowiregJ. We use index n-n'-n"-n"'u^ to describe the 
nanowire consisting of coaxial shells with n, n', n", n'" 
helical atomic rows [n > n' > n" > n'") each. Thus, 
the nanowire structures in Fig.l can be defined as 8-2, 
9-3, 9-4, 5-1, 6-1, 12-6-1, and 17-12-6-1 respectively. The 
structure of the thinnest 5-1 wire is centered pentagonal, 
with the outer shell containing five strands and the inner 
shell as a single atomic row. The 6-1, 12-6-1, and 17-12-6- 
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1 wires in Fig.l constitute a hexagonal growth sequence 
containing two-shell, three-shell, and four-shell with a 
central single atomic row. In stead of central atomic 
row based structures, the 8-2, 9-3, 9-4 wires are com- 
posed of eight or nine strands in the outer shell and two 
to four strands in the inner shell. Moreover, multi-shell 
structures such as 13-8-2, 14-9-3, 15-9-4 are obtained in 
our simulation. Our systematical illustration of various 
multi-shell growlh sequence of the metal nanowire can be 
found elsewherellj. 
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FIG. 1. Morphology of Ti nanowires with diameter from 
0.75 to 1.7 nm. Top view (left) and side view (right) are 
presented. These multi-shell structures can be identified as 
8-2, 9-3, 9-4, 5-1, 6-1, 12-6-1, and 17-12-6-1 respectively. 



well-reconstructed helical surface structures, the melting 
temperatures of nanowires are typically higher than those 
of clusters with comparable size (see Table I). 



1 

.4-1 
CS 

> 



U 

a 



S 
o 



U 

Qi 

a 



-3.80 
-3.85 
-3.90 
-3.95 
-4.00 
-4.05 

0.35 
0.30 
0.25 
0.20 
0.15 
0.10 
0.05 
0.00 



-3.85 
-3.90 
-3.95 
-4.00 
-4.05 
-4.10 
-4.15 
0.35 
0.30 
0.25 
0.20 
0.15 
0.10 
0.05 
0.00 



T ' 1 ' 1 ' 1 ' 1 ' r— 



—\ ' 1 ' 1 ' r- 



400 600 800 1000 1200 1400 

Temperature (K) 



400 600 800 1000 1200 1400 

Temperature (K) 



Starting from those optimized structures at low tem- 
perature, we simulate the melting behavior of the tita- 
nium nanowires using constant-temperature MD method. 
As representative examples. Fig. 2 (a), (b), (c), (d), (e) 
present the temperature dependence of the internal en- 
ergy E, mean square bond length fluctuation S, and ther- 
mal capacity Cy for 5-1, 6-1, 9-4, 12-6-1, and 17-12-6-1 
nanowires respectively. From those curves, we can define 
the overall melting temperature T™ of the system as the 
temperature at half of saturated S in liquid states'. For 
the purpose of comparison, we have also simulated the 
melting of some titanium nanoclusters with equilibrium 
structures such as icosahedral TissEH. Table 1 summa- 
rizes the Tm of various nanowires with different multi- 
shell structural pattern from our simulations. In general, 
the melting temperatures of titanium nanowires are much 
lower than the bulk value (1943 K) and increase with the 
nanowire size. Because of the 1-D periodicity and the 
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FIG. 2. The rms bond length fluctuation 5, internal energy 
and heat capacity C„ as functions of nanowire temperature 
for various nanowires. (a): 5-1; (b): 6-1; (c) 9-4; (d) 12-6-1; 
(e) 17-12-6-1. 

Table L Overall melting temperatures T^ for titanium 
nanowires and clusters. The diameter D of the nanowires 
or cluster are given. 
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To explore the size dependence of nanowire melting 
temperature, we plot the overall melting temperature Tm 
versus wire diameter (1/D) in Fig. 3. It is found that the 
T,„ for the 6-1, 12-6-1, 17-12-6-1 wires under the same 
hexagonal growth sequence fit well to a linear dependence 
as: 

where Tq = 1542 K is the extrapolated infinite limit, 
r] = 682 K • nm describes the linear dependence of Tm 
with diameter D (in unit of nm). Such 1/D dependence 
of melting temperature for nanowire is similar to-the-|Well- 
known size relationship for metal nanoclustersc3~E3. As 
compared with the linear dependence in the hexagonal 
wires (6-1, 12-6-1, 17-12-6-1), 8-2 and 13-8-2 wires al- 
most belong to the same size relationship while the Tm 
of other nanowires deviate from such linear fit. These dif- 
ferences indicate that the atomic structures of nanowires 
play significant role in determining the melting behavior 
of nanowires. Among those wires, the T„ of 5-1 wire is 
much higher than the 6-1 wires with comparable size. On 
the other hand, the nanowires with three or four atomic 
strands in the internal shell like 9-3, 14-9-3, 9-4, 15-9- 
4, have lower melting temperature than the nanowires 
with one atomic strand in the center, e.g., 12-6-1. This 
effect might be understood by the relatively looser inter- 
nal structures of these 9-3 or 9-4 based nanowires (see 
Fig.l). Therefore, we suggest that the multi-shell metal 
nanowires with tightly internal structures should have 
higher thermal stability. 
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FIG. 3. Overall melting temperature Tm (K) vs. nanowire 
diameters 1/D (1/nm) for various titanium nanowires with 
multi-shell structures. Dashed line is linear flt for hexagonal 
wires, i.e., 6-1, 12-6-1, 17-12-6-1. 

We now discuss the nanowire melting behavior by ex- 
amining the temperature dependent curve of internal en- 
ergy E, bond fluctuation {5) and thermal capacity {Cy) 
with different size and geometries. As shown in Fig.2(a- 
e), the melting of 5-1, 6-1 and 9-4 wires all start from 
about 500 — 600 K. Their difference in the internal struc- 
tures are only reflected in the slightly different variations 
of rms bond length fluctuation 5 with temperature. Sim- 
ilar effect has been also found in other thinner nanowires 
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like 8-2, 9-3. For all the thinner wires (e.g., 5-1, 6-1, 
9-4 in Fig.(a-c)), we haven't seen clear peak in the C„ 
curve in Fig. 2, which is a characteristic for the first- 
order solid to liquid phase transition. Accordingly, the 
increase of internal energies is also relatively smooth in 
these small nanowires. Only the dramatic jump in bond 
length fluctuation d with increasing temperature demon- 
strates the occurrence o£jnelting. This behavior is similar 
to that in small clusterscil. It can be undeptood by finite 
size effect in first-order phase transitionEjO, which leads 
to the coexistence of solid-apd liquid states in a rather 
broad temperature regionB9TH3. On the other hand, the 
thicker wires, such as 12-6-1 and 17-12-6-1, show clear 
characteristics of the solid to liquid phase transition (see 
Fig.2(d,e)). In particular, there are sharp peaks in the 
thermal capacity curves at 1050 K and 1100 K for 12-6-1 
and 17-12-6-1 wires respectively. 
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FIG. 4. Snapshot of the 17-12-6-1 nanowire at different 
temperatures during the melting process. A structural trans- 
formation is clearly seen at 1050 K. 

To further understand the melting process in the tita- 
nium nanowires, we have traced the structural change of 
the 17-12-6-1 wire with the temperatures in the melting 
region. Fig. 4 gives the snapshot of the nanowire at differ- 
ent temperatures. One can see that the titanium atoms 
only weakly oscillate at their equilibrium positions be- 
fore the overall melting (950 — 1000 K). Although there 
are small deformation on the sharp of nanowire and the 
diffusion of interior atoms, the top view of the nanowire 
cross-section still shows hexagonal structure with distinct 
atomic shells. It is interesting to find a structural trans- 
formation of the nanowire at 1050 K. In additional to 
the thermal enhanced 1-D diffusion along wire axis, the 
atomic arrangements of the nanowire have transformed 
from cylindrical helical to bulk-like rectangular. In the 
further melting process, the structures of nanowire re- 
mains bulk-like. The structural transition from cylindri- 
cal helical to bulk-like rectangular has also been observed 
in the size evolution of gold najipwire optimized struc- 
tures in our previous simulationta. 

From the above results and discussions, we make the 
following conclusions. (1) The thermal stability of tita- 
nium nanowires significantly dependent on the structures 
and sizes of the wire. (2) The melting temperature of ti- 
tanium nanowire is lower than the bulk's, but higher than 
titanium cluster's with comparable size. Linear depen- 



dence of melting temperature upon inverse of nanowire 
diameter {1/D) is obtained. (3) The coexistence mecha- 
nism of the solid state and the liquid state due to finite 
size effect plays a significant role in the melting process of 
the thinner nanowires {D < 1.2nm). (4) For the thicker 
nanowires, i.e., D ^ 1.7 nm, we find the structural trans- 
formation from the helical multi-shell cylindrical to the 
bulk-like during the melting process. 
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